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Abstract — This research formulates some dispersion equations with formulating procedure similar to the 
one in formulating dispersion equation of the small amplitude and long wave theory, i.e. by applying 
velocity potential equation on the Bernoulli surface equation and Kinematic Free Surface Boundary 
Condition equation. 

Furthermore, this research uses non-linear term of the Bernoulli equation, whereas the Kinematic Free 
Surface Boundary Condition equation is applied with two scenarios, i.e. neglected non-linear term and not- 
neglected non-linear term 

Wave length from various dispersion equations that are obtained are then compared with breaker length of 
the breaker index equation. 

This research aims only to show that using similar governing equations can be obtained some dispersions 


equations to produce different wave length. 

Keywords — Dispersion Equation, Wave Length. 

I. INTRODUCTION 

Wave length is an important parameter of a water wave. 
Various phenomena in a water wave that are determined by 
wave length are among others shoaling and breaking, 
refraction and diffraction, wave force on a structure, and 
sediment transportation by a wave. Therefore, a dispersion 
equation that produces an appropriate wave length is 
needed. 

Dean (1991) formulated dispersion equation an equation to 
calculate wave length, using two basic equations, i.e. 
Bernoulli equation at the surface and Kinematic Free 
Surface Boundary Condition (KFSBC) equation. In both 
equations, the non-linear term is neglected, by applying an 
assumption of small amplitude and long wave. 

In this research, dispersion equation is formulated using 
similar governing equation, i.e. Bernoulli surface equation 
and KFSBC equation by keep applying the small amplitude 
wave assumption but without applying the long wave 
assumption. The non-linear term at the Bernoulli equation 
is still used, whereas the KFSBC is applied in two 
scenarios, i.e. the neglected non-linear term and the not 
neglected non-linear term. In the scenario where KFSBC 
equation is not neglected, the formulation is applied with 
two different approaches. 


The wave length from the resulting dispersion equation is 
compared with breaker length calculated by breaker index 
equations from Komar and Gaughan (1972), Me. Cowan 
(1894) and Miche (1944). Breaker height is calculated 
using equation from Komar and Gaughan (1972), using 
input breaker height breaker depth is calculated 
withMc.Cowan (1894) equation, using input breaker height 
and breaker depth, breaker length is calculated using Miche 
(1944) equation. 

II. VELOCITY POTENTIAL EQUATION 

By completing Laplace equation with variable separation 
method. Dean (1991) obtained the following velocity 
potential equation, 

<P{x,z,t ) = Acoskx ( Ce kz + De~ kz )sin(<rt) 
+Bsinkx(Ce kz + De~ kz )sin(ot) ....(1) 

This equation has two components, i.e. cos/cxcomponent 
and sinkx component. Hutahaean (2019) shows that both 
components have similar wave constant, where (1) can be 
written as, 

<P(x,z,t) — A (coskx + sinkx ) ( Ce kz + De~ kz ) 
sin(at ) (2) 
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Thecoskxandsinkxfu notions have an intersection point 
where the two functions have similar values, so at that 
point the velocity potential equation can be written as, 

<P(x,z,t ) = 2Acoskx ( Ce kz + De~ kz )sin(ot ) 

A new constant is defined, i.e. A — 2A 

0(x,z,t ) = Acoskx ( Ce kz + De~ kz )sin(ot) 

....(3) 

This velocity potential (3) is an equation at the 
characteristic point i.e. Appoint wherecosfcx = sinkx. The 
formulation of A, C,D constants using (3) will produce a 
constant value at the characteristic point, where the 
constant value applies at all points at the wave curve. 

At the formulation A, C, ^constants with (3), flat bottomis 
used (Dean (1991), i.e.by applying the bottom water 
kinematic boundary condition, where at the flat bottom 
^ = 0 applies where h(x) is water depth. Hence, the 
bottom water kinematic boundary condition in the form of 

w = — u— becomesw = 0, or— = Oat z — —h. The 

dx a z 

following is obtained 

0(x, z, t ) = Gcoskx coshQi + z)sinot .(4) 

Where wave constant G — 2ADe kh is defined. The detail 
of the formulation can be seen in Dean (1991). 

III. DISPERSION EQUATION OF THE LINEAR 
WAVE THEORY 

3.1. Water Surface Equation of the Linear Wave Theory 
In this section dispersion equation of linear wave theory or 
small amplitude and long wave theory is formulated with a 
procedure corresponds to the one in Dean (1991). 


30 

~~3t~ +9T1 = 

Substitute the potential flow equation, water surface 
equation is obtained, i.e.. 


Go C(t) 

n(x,t ) = —coskx cosk(h + ri)cosot H- 

9 9 

or 


Go / m 

ri(x,i) = —coskx coskh 1 + —) cosot + 
a V hs 


m 


For a very small wave amplitude A, then^ « 1, therefore 

the last equation becomes. 

Go C(t) 

ri(x, t) = —coskx coskhcosot H- 

9 9 

r\(x, t)has an average value against time that is equal to 
zero, then C(t) = 0, the water surface equation becomes, 
Gocoskh 

rj(x, t) =- coskx cosot 

9 


r]{x, t) = Acoskx cosot 
It is defined that 

Gacoskh 


.( 6 ) 


A = 


■(7) 


where 4is the wave amplitude. (6) can be written to be an 
equation for G, i.e. 


G = 


Ag 


acoshkh 


.( 8 ) 


0(x,z,t ) = Gcoskx coshk(h + z)sinot. 


■(9) 


(9) is velocity potential equation of the linear wave theory. 
The velocity particle equation in the horizontal direction at 
the surface for small amplitude and long wave is 


3.1. Applying Bernoulli Equation at the Surface, 

The Bernoulli equation at the surface is, 

- ^ + \ K 2 , + <) + grj + J = c(t).(5) 

By applying an assumption that the wave amplitude is very 
small then r l ;, where^is the water surface elevation vis-a-vis 
still water level, will also be very small, so that the surface 
pressure can be considered as equal on the entire surface 
and if the reference is the atmospheric pressure then p Tj — 
0. 

By applying an assumption that the wave amplitude is very 
small, then the velocity of particles uandw are also very 
small number so that in Bernoulli equation the second term 
is much smaller than the third term and therefore can be 
neglected. 


a<2n 

=-= Gksinkx coshkhsinot 

" z=v 

.( 10 ) 

Whereas the particle velocity in the vertical direction at the 
surface is 


a<2n 

w„ =-= —Gkcoskx sinhkhsinot 

v az z=v 

.(ii) 

To ease the writing, constant Gis still used with the value 
as in (8). 


3.2. Applying the Surface Kinematic Boundary Condition 

The next constant to be formulated its equation is wave 
number k, which will be formulated using surface 
kinematic boundary condition equation, i.e.. 
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a?? . a?; 

W„ =-h U„ — 

71 at v ax 


....( 12 ) 


For small amplitude andlong wave. 


—is a very small 
ax J 


number and can be neglected, hence (9) becomes 


w, 


...(13) 


— 

71 at 
Substitute (11) and (6), 


—Gkcoskx sinhkhsinot = — - 


Go 2 coskh 


9 


-coskx sinot 


A relation is obtained, i.e. 
er 2 = gk tanhkh ....(14) 

This equation is a dispersion equation of linear wave 
theory to calculate wave number k, whereas wave length 
can be calculated with a relation L 

k 

In deep water, the relation tankh — 1 applies, hence, 
dispersion equation becomes, 
a 2 = gk 


From this equation, wave number in the deep water is 
obtained, i.e. 


k 0 


a 


....(15) 


- k 2 osin 2 kx cosh 2 khsinotcosot 

9 
G 2 

- k 2 ocos 2 kx sinh 2 khsinotcosot 

9 

(17) 


4.2. Applying KFSBC Equation 

As with the formulation of dispersion equation in the 
previous section, dispersion equation is formulated using 
KFSBC equation, i.e. (12). Substitute (9) and (10) to (12), 

— = —Gkcoskx sinhkh — Gksinkx coshkh — ....(18) 
at ax y ' 

By equalizing (17) with (18) and applying it at the 
characteristic point, wher ecoskx = sinkx and the 
characteristic point is also applied at domain time t, 
i .e.cosot — sinot , the following is obtained 

(j2 G 

- coshkh - k 2 o ( cosh 2 kh + sinh 2 kh ) = 

9 2 g 

dp 

—k sinhkh — k coshkh — 

3x 

.(19) 


IV. DISPERSION EQUATION OF SHORT 
WAVE 

In this section, the formulation of dispersion equation will 
be done using Bernoulli equation at the surface and 
KFSBC equation as in the formulation of dispersion 
equation of linear wave theory. However, the long wave 
assumption is not applied and a complete Bernoulli 
equation is used instead. 

4.1. Water Surface Equation of a Complete Bernoulli 
Equation 

In this section, water surface equation will be formulated 
using a complete Bernoulli equation i.e. (5). Substitute (9), 
(10) and (11) to (5), usingC(t) = 0, 

Go 

n(x,t) — —coskx coshkhcosot 
9 

G 2 

- k 2 sin 2 kx cosh 2 khsin 2 ot 

2 9 

G 2 

- k 2 cos 2 kx sinh 2 khsin 2 ot ...(16) 


4.2.1. KFSBC linearandG linear 

In this section, the long wave assumption is applied where 
the second term in the right side (19) is very small 
compared to the first term, and relationGi.e. (8)isused to 
obtain dispersion equation in other form, i.e. , 

o 2 = gk tanhkh — ^-k 2 ( 1 + tanh 2 kh ) ...(20) 

Equation (20) shows that wave number or wave length is 
also determined by wave amplitude. 

4.2.2. KFSBC nonlinear, G linear, — linear 

ax 

In (19),Gof (8) is substituted to (19) whereas |pis 
substituted with (6) and was applied at the characteristic 
point to obtain 

o 2 — gk tanhkh — ^k 2 ( 2 + tanh 2 kh ) ...(21) 

(21) is another form of dispersion equation, where the 
wave amplitude is the parameter, so that the resulting wave 
length is also determined by wave amplitude. 


(16) is non-linear water surface equation for small 
amplitude wave. Time differential from (16) is, 


3 g Go 2 

— =- coskx coshkhsinot 

at g 


4.2.3. KFSBC nonlinear, Gand— nonlinear. 

SubstituteJpto (19) with differential (16) against 
horizontal x axisand apply to the characteristic point , the 
following equation is obtained 
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Gk 2 o ( cosh 2 kh + sinh 2 kh ) = 

—2o 2 coshkh + 2gk sinhkh 

— (Gok 2 cosh 2 kh + k 4 coskhj .(22) 

The water surface elevation r /of (6) at the characteristic 
point is. 



Then (16) is applied to the characteristic point and is 
equalized with (23), and the following equation is 
obtained, 

G 2 

gA — Go coshkh —— k 2 (cosh 2 kh + sinh 2 kh ) 

4 

.(24) 

(22) and (24) are two simultaneous equations with 
unknown wave constant Gand wave number k. With input 
wave amplitude A, water depth /land wave period 

2n 

T where <7 = —, wave contant Gand wave number /ccan be 

T 

calculated with (22) and (24). 

The result of wave length calculation with (14), (20), (21) 
and with the system of equation (22) and (24) is presented 
in table (1). On that table,L 14 is the wave length calculated 
with (14), L 2 q is the wave length calculated with (20), and 
so forth 

Table 1. Wave Length from four dispersion equations 
Wave length L (m) 


h 

(m) 

L 14 

(*20 

(*21 

(*22+24 

20 

88,79 

83,66 

80,43 

80,42 

19 

87,63 

82,57 

79,36 

79,39 

18 

86,35 

81,38 

78,19 

78,26 

17 

84,96 

80,07 

76,88 

77,01 

16 

83,45 

78,63 

75,45 

75,63 

15 

81,79 

77,04 

73,87 

74,11 

14 

79,98 

75,31 

72,12 

72,42 

13 

78,01 

73,4 

70,19 

70,56 

12 

75,85 

71,3 

68,07 

68,49 

11 

73,49 

68,98 

65,71 

66,21 

10 

70,9 

66,43 

63,1 

63,66 


Table (1), presents the result of wave length calculation 
with (14), (20), (21) and ((22)+(24)) using a wave with 
wave period T — 8seconds and wave amplitude A — 1.0 
m. The result of the calculation shows thatL 14 is the 
longest, whereas, L 21 andL 22+24 is more or less equal 
although L 21 is relatively shorter. In addition, there is a 
constraint at ((22)+(24)), i.e. it cannot be used in a shallow 
water. Henceforth, ((22)+(24)) can no longer be used. 

Wave length calculation is then done with wave period 
T — 8seconds and wave amplitude A = 1.0 m in a shallow 
water, with the result of the calculation as presented in 

table (2), whereG = Lli ~ L21 x 100%. 

l . 14 

Table 2. Comparison of wave length in a shallow water. 


h 

Wave 

length L 

(m) 

S 

(m) 

(*14 

(*20 

(*21 

(%) 

20 

88,79 

83,66 

80,43 

9,42 

18 

86,35 

81,38 

78,19 

9,46 

16 

83,45 

78,63 

75,45 

9,58 

14 

79,98 

75,31 

72,12 

9,83 

12 

75,85 

71,3 

68,07 

10,26 

10 

70,9 

66,43 

63,1 

11 

8 

64,9 

60,47 

56,94 

12,27 

6 

57,5 

52,98 

49,08 

14,65 

4 

48,01 

43,14 

38,38 

20,04 

2 

34,69 

28,5 

19,34 

44,24 


In the deep water, the difference between L 14 andL 21 could 
reach 9.5 %, whereas in shallow water the difference could 
reach 44 %. 

Furthermore, the effect of wave amplitude A on (21) will 
be studied using a wave with wave period of T — 8 sec., 
with various wave amplitudes, i.e. 0.20 m, 0.60 m and 1.0 
m, with the result of the calculation as presented in 
table(3). It shows that the larger the wave amplitude the 
shorter the wave length. It can be concluded that wave 
amplitude is to shorten the wave length. 

Table.3: Wave Length from (21) at various wave amplitude 
values 

Wave Length L (m) 

h 4 = 0.2 4 = 0.6 4 = 1.0 

(m) (m) (m) (m) 

20 87,28 84,04 80,43 
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18 

84,88 

81,72 

78,19 

16 

82 

78,91 

75,45 

14 

78,56 

75,52 

72,12 

12 

74,45 

71,44 

68,07 

10 

69,5 

66,49 

63,1 

8 

63,48 

60,42 

56,94 

6 

56,02 

52,8 

49,08 

4 

46,37 

42,75 

38,38 

2 

32,6 

27,66 

19,34 

To see the effect of the difference in wave length, particle 
velocity in the direction of horizontal xis used with the 
result of the calculation as presented in table (4), using a 
wave with wave period T — 8 sec., wave amplitude A = 


1.0 m, and velocity calculated at z — —0.25 h. The 
calculation of the wave number is done using (14), (20) 
and (21). In table (4) t ul4is the particle velocity 
calculated using wave number from (14), and so forth, 
whereas# = u21 ul4 xl00 %. 

u 14 

Table. 4. Particle velocity in the direction of horizontal x, u 


h 

ul4 

u20 

u21 

5 

(m) 

(m/sec) 

(m/sec 

(m/sec 

(%) 

20 

0,66 

0,68 

0,69 

5,66 

18 

0,7 

0,72 

0,74 

6,02 

16 

0,74 

0,77 

0,79 

6,47 

14 

0,8 

0,83 

0,86 

7,05 

12 

0,87 

0,91 

0,94 

7,84 

10 

0,96 

1,01 

1,04 

8,97 

8 

1,08 

1,14 

1,19 

10,76 

6 

1,25 

1,34 

1,43 

13,98 

4 

1,55 

1,7 

1,88 

21,55 

2 

2,2 

2,64 

3,73 

69,64 


The difference betweeiml4 andu21is quite large where the 
shallower the water the greater the difference. 

V. COMPARISON WITH BREAKERINDEXES. 

As a comparator of wave length produced by dispersion 
equation, breaker length of breaker Indexesare used. The 
procedure of calculating the breaker length using breaker 
indexes is as follows 
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Breaker height that is calculated with the Komar and 
Gaughan equation (1972) is 

a = 0.56gi) _Vs .(24) 

n 0 \L,q/ 

H b is breaker height, H 0 is deep water wave height andL 0 is 
deep water wave length calculated with (15). 

Breaker depth that is calculated with McCowan (1894) 
equation is 

^ = 0.78 .(25) 

H 

h b is breaker depth, whereas breaker height Ilfs obtained 
from (24). 

Breaker length that is calculated using Miche (1944) 
equation is 

^ = 0.142tanh ....(26) 

L b is breaker length. Breaker height H b was obtained from 
(24) whereas breaker depth h b was obtained from (25). 

The calculation of breaker height with (24) requires an 
input of deep water wave height H 0 and wave period 7'for 
the calculation of deep water wave length L 0 . Those two 
parameters were obtained by applyingWiegel equation 
(1949,1964). 

By establishing a wave period T, deep water wave height 
H 0 is calculated using Wiegel equation (1949,1964), i.e. 

T = 15.6 P .(27) 

a] 9 
or 


Wheregis gravitational velocity ^9.81 m / sec 2 ^> 

deepwater wave height H 0 is in meter unit.The result of the 
calculation of the breaker length with this procedure is 
presented in table (5). 


Table.5: Breaker length L b in various wave periods. 


T 

Wo 

H b 

h b 

4 

(sec.) 

(m) 

(m) 

(m) 

(m) 

6 

1,45 

1,69 

2,16 

20,41 

7 

1,98 

2,3 

2,95 

27,78 

8 

2,58 

3 

3,85 

36,28 

9 

3,27 

3,8 

4,87 

45,92 

10 

4,03 

4,69 

6,01 

56,69 
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11 4,88 5,68 7,28 68,6 

12 5,8 6,75 8,66 81,63 

Furthermore, with breaker depth h b and with an assumption 
of sinusoidal wave where A = -y-, wave length is calculated 
with (14), (20) and (21), with the result as presented in 
table (6), where on the table,L 14 is wave length of (14), 
L 20 is wave length of (20) andL 21 is wave length of(21). 

Table.6: The comparison of wavelength of dispersion 
equation with breaker length L b 


T 

4 

Ll4 

4o 

4i 

(sec.) 

(m) 

(m) 

(m) 

(m) 

6 

20,41 

26,53 

22,05 

15,67 

7 

27,78 

36,11 

30,02 

21,33 

8 

36,28 

47,17 

39,21 

27,86 

9 

45,92 

59,7 

49,62 

35,26 

10 

56,69 

73,7 

61,26 

43,54 

11 

68,6 

89,18 

74,13 

52,68 

12 

81,63 

106,13 

88,22 

62,7 


Table (6) shows that L 20 is the closest to L b . However, if it 
is viewed based on wavesteepness criteria of Michell 
(1893) where^ = 0.142, then the one that makes wave 
steepness to be closer to critical wave steepness isL 21 , as 
presented in table (7). 


Table. 7: Comparison of wave steepness 


T 

H b 



Hb_ 

(sec.) 

4 

L 14 

4o 

4i 

6 

0,083 

0,064 

0,077 

0,108 

7 

0,083 

0,064 

0,077 

0,108 

8 

0,083 

0,064 

0,077 

0,108 

9 

0,083 

0,064 

0,077 

0,108 

10 

0,083 

0,064 

0,077 

0,108 

11 

0,083 

0,064 

0,077 

0,108 

12 

0,083 

0,064 

0,077 

0,108 


VI. 

CONCLUSION 



As a conclusion, from a governing equation can be 
obtained e some dispersion equations that produce various 
wave lengths. The higher the level of the precision, the 
shorter the wave length. Even using an assumption of small 
amplitude wave, dispersion equation with wave amplitude 


as its parameter can be resulted. The influence of wave 
amplitude is to shorten the wave length. 

Variety of dispersion equations producing variety of wave 
lengths require a criteria on the appropriate wave length. 
One of the criteria that can be used is critical wave 
steepness. 

Further research needed is formulating dispersion without 
applying an assumption of small amplitude wave and by 
taking into account the criteria of critical wave steepness. 
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